{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "favorite-pressing",
   "metadata": {},
   "source": [
    "# Solving Transcendental Equations using fsolve()\n",
    "*March 19, 2021*\n",
    "\n",
    "- First import some required modules"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "outside-kernel",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from scipy.optimize import fsolve\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "prime-punch",
   "metadata": {},
   "source": [
    "The transcendental equation that we will solve in this example is:\n",
    "\n",
    "$\\ln\\left(\\dfrac{\\dot{Q}RT}{\\dot{V}P_0L_\\mathrm{m}}\\right) = \\dfrac{L_\\mathrm{m}}{R}\\left(\\dfrac{1}{T_0}-\\dfrac{1}{T}\\right)$\n",
    "\n",
    "This equation determine the equilibrium temperature of a boiling liquid that is being cooled by evaporative cooling.  It's vapor pressure is reduced via a pump with a pumping speed $\\dot{v}$ and the liquid is absorbing a heat load of $\\dot{Q}$ from its surroundings.\n",
    "\n",
    "The dependence of the boiling temperature on the vapor pressure is determined from the Calusius-Clapeyron relation.  The liquid has an equilibrium temperature of $T_0$ when the vapor pressure is $P_0$.\n",
    "> $R$ is the universal gas constant<br>\n",
    "$L_\\mathrm{m}$ is the molar latent heat\n",
    "\n",
    "- Enter some numerical values of the various quantities.  The goal is to find the equilibrium temperature which cannot be solved for analytically."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "sitting-climate",
   "metadata": {},
   "outputs": [],
   "source": [
    "Qdot = 2.5 # W\n",
    "vdot = 1e-3 # m^3/s\n",
    "Lm = 90 # J/mol\n",
    "R = 8.31 #J/K mol\n",
    "P0 = 99.23e3 # Pa\n",
    "T0 = 4.2 # K"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "after-envelope",
   "metadata": {},
   "source": [
    "- Create a function which defines our transcendental equaition.  We want to find $T$ that makes \"Fcn(T)\" equal to zero."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "sticky-plaintiff",
   "metadata": {},
   "outputs": [],
   "source": [
    "def Fcn(T):\n",
    "    return np.log(Qdot*R*T/(vdot*P0*Lm)) - (Lm/R)*(1/T0 - 1/T)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "exact-contribution",
   "metadata": {},
   "source": [
    " - We use *fsolve()* to find $T$.  *fsolve()* requires two arguments.  The first is the function that we defined and the second is an initial guess at the solution.  If your guess is too far from the actual solution, *fsolve()* may return an error."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "graduate-johnston",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([1.29141355])"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sol = fsolve(Fcn, 2)\n",
    "sol"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "still-browse",
   "metadata": {},
   "source": [
    "The equilibrium temperature for a heat load of $\\dot{Q}=2.5~\\mathrm{W}$ and a pumping speed of $\\dot{v}=1.0\\times 10^{-3}~\\mathrm{m}^2/\\mathrm{s}$ is $T=1.291~\\mathrm{K}$.\n",
    "\n",
    "We can look to see if this solution makes sense.  We do this my plotting the left-hand and right side of the transcendental equation as a function of $T$ on the same graph and look for the intersection.\n",
    "\n",
    "- First define functions for the LHS & RHS of the transcendental equation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "marked-doubt",
   "metadata": {},
   "outputs": [],
   "source": [
    "def LHS(T):\n",
    "    return np.log(Qdot*R*T/(vdot*P0*Lm))\n",
    "\n",
    "def RHS(T):\n",
    "    return (Lm/R)*(1/T0 - 1/T)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "excessive-bristol",
   "metadata": {},
   "source": [
    "- Next, make a list of some temperature values near the solution that *fsolve()* found. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "minimal-translation",
   "metadata": {},
   "outputs": [],
   "source": [
    "TT = np.arange(0.9, 1.7, 0.01)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "honest-thesaurus",
   "metadata": {},
   "source": [
    "- Now, plot the two curves and visually inspect the intersection."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "elegant-enclosure",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD4CAYAAADxeG0DAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAf80lEQVR4nO3deXxU9bnH8c8DiIAgO7KJoIhLqyJGq1TcoNaKQFErouJe3OBavW4V11q1F/cqihRT96LWFVEURKEFqQYQpWDF1uUCsq+iEEJ+948nczPEhCTMZM6ZzPf9ep3XnDMzzHkY4JvD7/wWCyEgIiLZq07UBYiISGoU5CIiWU5BLiKS5RTkIiJZTkEuIpLl6kVx0latWoXOnTtHcWoRkaw1a9aslSGE1mWfjyTIO3fuTEFBQRSnFhHJWmb2VXnPq2lFRCTLKchFRLKcglxEJMspyEVEspyCXEQkyynIRUSynIJcRCTLKchFRGpYCDB7Nlx9NSxfnv7Pj2RAkIhILvjXv+Avf/Hts8+gXj046ijo1y+951GQi4ik0eLFMG4cPPusX4WbwdFHw3//N5xyCrRsmf5zKshFRFK0di28+CI88wy89543peTlwT33wKBB0KFDzZ4/LUFuZlcBdwGtQwgr0/GZIiJxtnkzvPGGh/frr/tx165w001wxhnQrVvmakk5yM1sd+BnwNeplyMiEl/FxTBjBjz1FDz/vF+Jt2kDF10EZ54Jhx7qTSmZlo4r8vuAa4BX0/BZIiKxs3Chh/dTT8GXX0KjRnDyyXDWWdC7t9/EjFJKpzez/sDiEMJcq+THkJkNBYYCdOrUKZXTiojUuNWr4bnn4MknYeZMqFPHQ/t3v4OBA6Fx46grLFVpkJvZZKBtOS+NAK4Hjq/KiUIIY4AxAHl5eaEaNYqIZMSWLTBxIjzxBIwfD4WF8OMfw8iR3u5d0zctd1SlQR5C6FPe82Z2ANAFSFyNdwRmm9lhIYSlaa1SRKQGffwxPP6437hcvhxat4ZLL4Wzz4bu3aNp966OHW5aCSF8ArRJHJvZl0Ceeq2ISDZYtcqD+/HHYc4c2Gkn6N8fzjkHTjjBj7OF+pGLSM4oKoK334b8fHjtNW9K6dEDHnwQBg+umcE6mZC2IA8hdE7XZ4mIpNPChfDnP3vb95Il3nQybBicey4ceGDU1aVOV+QiUitt3Ah//Ss89hj87W/e6+TEE+Ghh6BvX6hfP+oK00dBLiK1RggwaxaMHetznWzYAHvvDXfe6Tcu27ePusKaoSAXkay3Zo3fuBw7FubOhYYN4Ve/ggsugF694t/rJFUKchHJSiH4cPkxY3y4/KZNcPDB8PDD3ue7adOoK8wcBbmIZJU1a3y05ZgxMH8+NGniXQZ//Ws45JCoq4uGglxEYi8EeP99ePTR0qvvww7zppRBg+I1XD4KCnIRia316+Hpp2H0aPjkE7/6Pu88GDrUR1yKU5CLSOzMmePh/cwz3o2wRw/405/g9NN19V0eBbmIxMKmTd5s8sgjPttgw4Y+2vKSS3y1HamYglxEIvXFFx7e+fk+/0m3bnDffX4Ds3nzqKvLDgpyEcm44mKf82TUKJgwwUdd9u8Pl10Gxx1X+/t9p5uCXEQyZt06n21w1Cif/6RNGxgxwpdK69gx6uqyl4JcRGrcggU+w+CTT/rNyyOOgFtugVNOgZ13jrq67KcgF5EaUVzsq8z/8Y8waZIH9uDBPutgrg7cqSkKchFJq/Xr/cblQw/Bv//ty6PdfruPvGzdOurqaicFuYikxb//7c0n+fk+62DPnnDHHb5QcTattpONFOQissNCgGnTvLvga69B3bo+ZP7yy+HQQ6OuLncoyEWk2goL4bnnPMDnzPEl0q6/3hcsrq1zfseZglxEqmz1ap918MEHfcm0/fbziayGDPGRmBINBbmIVOo///Gr7/x8+O476NPHZx484QQN3okDBbmIVOiDD+Duu+HFF739e/BguPJKOOigqCuTZApyEdlGov/3XXf5jcymTeHqq2H4cO9KKPGjIBcRwG9gPvusB/j8+bD77nDvvXDhhT4PuMSXglwkx23Y4Dcw77sPFi+GAw/0xRxOO039v7OFglwkRy1f7sPnR42CtWvh2GPhscfg+ON1AzPbKMhFcsxXX/kNzLFjYfNmH3l57bW+BqZkpzqp/GIzu8XMFpvZRyXbiekqTETSa8ECX6yha1fv+33GGf7ciy8qxLNdOq7I7wsh3J2GzxGRGjB7ts958tJLPmhn2DDvQrj77lFXJumiphWRWmr6dPj972HiRO9COGKEz4HSqlXUlUm6pdS0UmKYmX1sZvlmVuEKe2Y21MwKzKxgxYoVaTitiJQVAkyZ4jcujzwSCgr8avyrr+C22xTitZWFELb/BrPJQNtyXhoBzARWAgG4DWgXQji/spPm5eWFgoKC6lcrIuUKwdfA/N3vYMYMaNcOrrnG5wDfZZeoq5N0MbNZIYS8ss9X2rQSQuhTxRP8CXh9B2oTkR0UArz5pgf4P/7h7d6jRsH550ODBlFXJ5mSaq+VdkmHA4F5qZUjIlURArz+uvc26dsXli3zQT2ff+5TySrEc0uqNztHmll3vGnlS+CiVAsSkYqFABMm+MLFs2ZBly4+iGfIEI3CzGUpBXkIYUi6ChGRioXgE1ndcovfwNxzT59S9qyzFOCSnl4rIlJDQoC33oLDD4eTToKVK/0K/NNP4bzzFOLiFOQiMfXuu9Crly/esHSpt4F/9pnfyFSASzIFuUjMzJgBvXvDccfBl196L5TPPvOuhApwKY+CXCQm5szxHig//SnMmwf331/aC2XnnaOuTuJMQS4SsX/9CwYNgh494P334c47fY3Myy9XN0KpGs21IhKRr7+GW2+Fxx/3yaxuvNEns2rWLOrKJNsoyEUybOVKv+oeNcp7pfzXf8Fvfwtt2kRdmWQrBblIhnz7rS+ndtddsHGjzw1+yy3QqVPUlUm2U5CL1LAtW3w1nltv9aH0AwfC7bfDfvtFXZnUFgpykRoSAvz1r3D99d77pFcveOUVH9wjkk7qtSJSA6ZNgyOO8JXod97ZJ7iaOlUhLjVDQS6SRgsWwIABcPTR8L//68Pp5871/uFamV5qippWRNJg+XK/cTlmDDRq5G3gv/mN70sOKSyEDRtg/fptt+TnTj0VOndO62kV5CIp+P5774nyhz/Ad9/BxRfDzTdD69ZRVybVUlgI69aVhm3yftmtotfWrYPNmys/1/77K8hF4qC4GMaNg+uu8yaU/v1h5EjYZ5+oK8sxxcXerzMRruvWlW7Jobu9x6oGcL16vop1kyb+uOuuvqZet26lx8mvlXe86641svaeglykmmbMgCuugA8+gIMPhiefhGOOibqqLFRc7E0Oa9duG8BV2ZLDuJJ1h6lTpzREE6G6224ewGWfTw7qsuHcoEFsb3QoyEWq6Kuv/Ap83Dho396H1g8Z4jmRc0KATZs8TBNBvHbttvvJj+Xtb9hQeQjvtFNpqCa2rl1/+FwihMs73mWX2AZwuijIRSrx7bfwP/8Dd9/txzfd5CvUZ/Xq9CF4A38ifNeuhTVrfhjKyVvZ5wsLt3+OunU9SJs1Kw3VvfYq3U9+vqItxlfBcaIgF6lACPDMM3DttbBkCZxxhs+REpsh9Vu2bBvAZR+3t792rf/67dl5Z2je3AO3WTPf79Jl2wBOfr1sQOfAlXBcKMhFyvHhhz6N7PvvQ14evPAC9OxZAyfatMkDNrGtXr1tICe25OPE/saN2//s+vVLA7hZM2jRwhf7TA7f5C05mBNXw5IVFOQiSZYv95kI8/P9flh+vk9utd128C1btg3h5P3k58q+tmaNB/n2NGlSGsbNm3v7cGI/+Uq57HGzZj43rq6Ic4KCXATYUhgYdX8hN/9+J77/Hq46bRE39vuIXTcth7tXe/iW3RKh/O232//wJk38ajgRuPvuW7q/va1ZM+/yJlIJ/S2R2iVxE2/VKg/Zih4T26pVTPlmP4av/R3zw/78nIk8wOXs8/xn8HzS59av72Gc2Dp1gu7dSwO67GNiX2EsGaC/YRJfRUWl4Vvelvxa8v72Bnc0bAgtW0Lz5izaZR+uWn8Hz605ki5NV/HKgJfof+RqrOWdpUGc2Bo1UjOFxJaCXDJj82YP2ZUry38sb3/duoo/r149D+TEttdecNhhvp8I38R+8nMNG1JYCA884PODb93qc6Rcc01LGjY8OWNfh0g6Kcil+oqKSsO2qtv22pEbN4ZWrbYN5eTj8rbGjXfoCnnqVF+Vfv586NfPV6rfc88d/ypE4iDlIDez4cAwoAiYEEK4JuWqJHNC8G5sK1b4tnJl6X7iOPFcYn/Nmoo/r0kTD+HEts8+PoNUIpjLPrZs6f2Va9jSpXD11fD00z5f0fjxcNJJNX5akYxIKcjN7FhgAHBgCGGzmWn52KiF4PNPrFjhfemSQ7nscSKcK+oCt9NOHsKJID7kEH9MHCfCODmoMxDK1bF1Kzz6qK/S8/33cMMN3r1Q08tKbZLqFfklwB9CCJsBQgjLUy9JfuC77zyEy9uSAzrxWNHQ6caNS4O5fXs46KDS40QYJ+83aZLVN/hmz/ZpZT/8EHr39lXrNTuh1EapBnk3oJeZ3Q5sAq4KIXyYelm1XHGxtzEvX+6r8S5bVhrMyfuJraIRfLvsAm3aePB26OBT8SWCOPF88tawYWZ/nxHZsAFuvBEefNB/Hj3zDAwenNU/k0S2q9IgN7PJQNtyXhpR8uubA4cDhwLPm9meIfxwSjMzGwoMBegUm8kq0mjrVr8aTgTzsmXeMJt8nAjpFSv8/WXVreuBu9tu/ti1a+n+brt5OCcHdFbP2lQzXn0Vhg2DxYvhoot8bpRmzaKuSqRmVRrkIYQ+Fb1mZpcAL5UE9wdmVgy0AlaU8zljgDEAeXl5lcxdGRPFxd4/eenS0lAuu594XLGi/Ck5GzYsDeE99vAuconjso8tWuTonKipW7QIhg/3VeoPOACef94XPxbJBak2rbwCHAe8Z2bdgPrAylSLqnGbNsE333gAf/NN6X7iODmsi4p++OsbNPDwbdvWu0AcfrgfJ55L3t/BbnJSNcXF8MgjfgOzqMiXXLvySr9PK5IrUg3yfCDfzOYBhcA55TWrZEQI3ji6ZElpOJe3LV3qs8eVVaeON1e0a+cBfMAB/pg4TgR027Y+ab3COXL//Cf8+tc+Q2GfPjB6tHdBF8k1KQV5CKEQOCtNtVTuk098W7KkNLCT98u7KdiggYdxu3bwox/5v/hEOCeeb9vWQ7xu3Yz9VmTHbd4Md9zh7d+77gpPPOEr9ehnq+Sq7BrZOXo0PPyw7zdq5F3o2rXz/s2JUG7fvvT5du18XmX9C681Zs6ECy7wkZlnnukr2GvFesl12RXk11zjd7Tat8/6Ps5SPRs3epfC+++Hjh1hwgQ48cSoqxKJh+wK8j32iLoCicC778KFF8J//uPzpCSaVETEqa+bxNb69XDJJXDccX4veupUH52pEBfZloJcYmnSJO849Oij3p1w7lw46qioqxKJJwW5xMr69T4i8/jjfSzV9Olwzz2a5EpkexTkEhvvvONX4WPH+pSzc+ZodKZIVSjIJXIbN8Jll3kX/wYN4O9/h5Ejc2aOL5GUKcglUtOn+2y6jzwCV1wBH32kq3CR6lKQSyQ2bfJhAb16+Xwp770H996rq3CRHZFd/cilVpg7F846C+bN8xubd9/tc4uJyI7RFblkzNatPjvhoYf6CnNvvOGzLijERVKjK3LJiC++gLPP9huZv/qVt4m3bBl1VSK1g67IpUaF4LMTHnQQfPwxPPUUPPecQlwknRTkUmNWrYLTToNzz/XlRD/+2NvGNdeZSHopyKVGTJniV+Gvvurt4lOmaM4zkZqiIJe0KiyE667zwT2NG/v84ddeqzU7RGqSbnZK2ixcCIMHw6xZMHSo9wvfZZeoqxKp/RTkkrIQ/CbmpZfCzjvDSy/BwIFRVyWSO9S0IilZv97XyzznHMjL88E+CnGRzFKQyw4rKIAePWDcOLjtNp+9sGPHqKsSyT0Kcqm2EHzR4549YcsWX7nnhht0Q1MkKmojl2pZtcr7hb/+Ovzyl/DYY9CiRdRVieQ2XZFLlU2fDt27w9tvwx//6Dc1FeIi0VOQS6WKi32hh6OP9l4pM2bA8OEaoSkSF2pake1atcp7pEyYAKee6suwNW0adVUikiylIDez54B9Sg6bAWtDCN1TrEli4h//8LlSli6FBx/05dh0FS4SPykFeQhhUGLfzO4B1qVckUQuBBg1Cq68Ejp08LbxvLyoqxKRiqSlacXMDDgNOC4dnyfR2bDBh9ePGwcnneRT0OqGpki8petmZy9gWQhhYZo+TyKwYAEcdhg8/zzceafPXKgQF4m/Sq/IzWwy0Lacl0aEEF4t2R8M/KWSzxkKDAXo1KlTNcuUmvbCC3D++dCoEUyeDMceG3VFIlJVFkJI7QPM6gGLgUNCCIuq8mvy8vJCQUFBSueV9NiyxaedvfdeOOIID/QOHaKuSkTKY2azQgg/uGOVjjbyPsCnVQ1xiY9ly2DQIB9iP3y4r2Zfv37UVYlIdaUjyE+nkmYViZ8PPoCTT4bVq30K2rPOiroiEdlRKQd5COHcNNQhGfTYYz53ePv2Pkqze/eoKxKRVGiIfg4pLIRLLoELL4SjjvJpaBXiItlPQZ4jli2D3r1h9Gi45hqYOBFatoy6KhFJB821kgMKCnzVnlWr4NlnfV1NEak9dEVeyz37LPTqBXXq+FB7hbhI7aMgr6WKi+G3v4Uzz/TRmh9+CAcfHHVVIlIT1LRSC23Y4AE+fjxcdJEvAqH+4SK1l4K8lvniC+jXDz79FB56yLsZaupZkdpNQV6L/O1vPshn61Z46y3vpSIitZ/ayGuJJ57w4G7RwheEUIiL5A4FeZYrLvZJr8491wf5zJwJe+8ddVUikklqWsli330HQ4b4avYXX+w3NXfaKeqqRCTTFORZaulS6N/fB/vcdx9cfrluaorkKgV5Fpo3D/r2hZUr4ZVXPNBFJHepjTzLTJoEPXtCUZH3UlGIi4iCPIvk58OJJ0KXLt4zpUePqCsSkThQkGeBEODGG+GCC+C44/xKvGPHqKsSkbhQG3nMFRZ6gD/9tD8+8oh6pojItnRFHmPr13tTytNPw223wZ/+pBAXkR/SFXlMLVkCv/gFzJ8Pjz8O55wTdUUiElcK8hiaP99DfPVqmDABjj8+6opEJM4U5DEzYwacdJJPOzt1qnqmiEjl1EYeI+PH+2RXrVrB++8rxEWkahTkMZGf7+tq/vjHviRbly5RVyQi2UJBHrEQ4M47vWth797w7rvQunXUVYlINlGQR6i4GK66Cq6/3hdFHj8eGjeOuioRyTa62RmRoiK48EJfEGLYMHjgAV/pXkSkuhQdEfj+ezjlFA/xW2/1ecQV4iKyo1K6Ijez7sBooAFQBFwaQvggDXXVWuvX+4yF06b54siXXRZ1RSKS7VJtWhkJ3BpCeNPMTiw5PiblqmqpVat8oM/s2fDUU3DmmVFXJCK1QapBHoBdS/abAktS/Lxa65tv4Gc/g88/h5dfhn79oq5IRGqLVIP8N8BbZnY33t7es6I3mtlQYChAp06dUjxtdvnyS+9auHw5vPkmHHts1BWJSG1SaZCb2WSgbTkvjQB6A1eEEF40s9OAx4A+5X1OCGEMMAYgLy8v7HDFWeazzzzEN26EyZPhJz+JuiIRqW0qDfIQQrnBDGBmTwKXlxy+AIxNU121wrx50KeP9xd/7z048MCoKxKR2ijVTm9LgKNL9o8DFqb4ebXG7NlwzDFQt673UFGIi0hNSbWN/NfAA2ZWD9hESRt4rps5E044AZo1g3fegb32iroiEanNUgryEMLfgUPSVEut8Pe/+6o+u+0GU6bA7rtHXZGI1HYaT5hG773nV+Lt2/tc4gpxEckEBXmaTJ7sV+J77OGB3r591BWJSK5QkKfB22/7qj577+0h3ra8zpoiIjVEQZ6it97yuVP23ddvbGoucRHJNAV5CiZOhAEDYL/9PMRbtYq6IhHJRQryHTRxIvzyl7D//t4+3rJl1BWJSK7SwhI74O23tw3xFi2irkhEcpmuyKtpyhRvTtl3X5g0SSEuItFTkFfD1KneO6VrVzWniEh8KMiraPp06NsXunTRjU0RiRcFeRV88IGv7NOxo4d4mzZRVyQiUkpBXomPPoKf/9z7h7/zjgb7iEj8KMi3Y/58X56tSRMP8Q4doq5IROSHFOQVWLjQV/apV89DvHPnqCsSESmf+pGX4+uvfWWfoiLvqbL33lFXJCJSMQV5GcuWeXPKunXw7rs+6EdEJM4U5EnWrPEbm4sW+ejNgw+OuiIRkcopyEt8+63PJ75gAYwfDz/9adQViYhUjYIc2LwZTj7Z+4u/8AIcf3zUFYmIVF3OB/nWrTBkiM+b8uc/e6CLiGSTnO5+GAJceqlfhd9zD5x7btQViYhUX04H+Q03wJgxcP31cOWVUVcjIrJjcjbI778f7rgDLroIfv/7qKsREdlxORnkzz4LV1wBp5wCo0aBWdQViYjsuJwL8kmTvC386KPh6aehbt2oKxIRSU1OBXlBgfdK2X9/ePVVaNAg6opERFKXUpCb2UFm9r6ZfWJm481s13QVlm6ff+4Dflq1gjffhKZNo65IRCQ9Ur0iHwtcF0I4AHgZuDr1ktJvxQo44QQoLoa33oJ27aKuSEQkfVIN8n2AaSX7k4BTUvy8tNu40dfZXLzYh9536xZ1RSIi6ZVqkM8D+pfs/wrYvaI3mtlQMysws4IVK1akeNqqKSqC00/3tvFx4+CIIzJyWhGRjKo0yM1sspnNK2cbAJwPXGZms4AmQGFFnxNCGBNCyAsh5LVu3Tp9v4MKzwfDhsHrr8NDD8GAATV+ShGRSFQ610oIoU8lbzkewMy6AX3TUVQ6jBwJjz4K110Hl1wSdTUiIjUn1V4rbUoe6wA3AKPTUVSqnnvOA/z00+H226OuRkSkZqXaRj7YzD4DPgWWAH9OvaTUTJ8O55wDRx7psxnWyame8iKSi1KaxjaE8ADwQJpqSdnChd4W3qkTvPKKBvyISG6oNderq1ZB374+b8obb0DLllFXJCKSGbViYYnCQp8A6+uvYcoU6No16opERDIn64M8BLj4Ypg6FZ55Bnr2jLoiEZHMyvqmlZEj/abmzTfDGWdEXY2ISOZldZC/9JJ3Mxw82INcRCQXZW2Qz5njiyYffjjk52txCBHJXVkZ5EuXQv/+3jNF3QxFJNdl3c3OTZtg4EBYvdoH/+y2W9QViYhEK6uCPAQYOhRmzoQXX4Tu3aOuSEQkelnVtHLXXfDUU3Dbbb5km4iIZFmQ77EHnHcejBgRdSUiIvGRVU0rgwb5JiIipbLqilxERH5IQS4ikuUU5CIiWU5BLiKS5RTkIiJZTkEuIpLlFOQiIllOQS4ikuUshJD5k5qtAL7awV/eCliZxnLSKa61xbUuiG9tca0L4ltbXOuC+NZW3br2CCG0LvtkJEGeCjMrCCHkRV1HeeJaW1zrgvjWFte6IL61xbUuiG9t6apLTSsiIllOQS4ikuWyMcjHRF3AdsS1trjWBfGtLa51QXxri2tdEN/a0lJX1rWRi4jItrLxilxERJIoyEVEslxsg9zMTjCzf5nZ52Z2XTmvNzezl83sYzP7wMx+nKG68s1suZnNq+B1M7M/ltT9sZn1iEld+5rZ+2a22cyuykRN1ajtzJLv6mMzm2FmB8WkrgElNX1kZgVmdmQm6qpKbUnvO9TMtprZqXGoy8yOMbN1Jd/ZR2Z2UybqqkptSfV9ZGb/NLOpcanNzK5O+s7mlfyZtqjyCUIIsduAusC/gT2B+sBcYP8y77kLuLlkf1/gnQzVdhTQA5hXwesnAm8CBhwO/CMmdbUBDgVuB67K8J9nZbX1BJqX7P8iRt9ZY0rvIx0IfBqX76zkPXWBKcAbwKlxqAs4Bng9k3+/qlFbM2A+0KnkuE1caivz3n7AlOp8flyvyA8DPg8h/CeEUAiMAwaUec/+wDsAIYRPgc5mtltNFxZCmAas3s5bBgBPBjcTaGZm7aKuK4SwPITwIbClpmsp59yV1TYjhLCm5HAm0DEmdX0bSv5lAbsAGesZUIW/ZwDDgReB5TVfkatiXZGoQm1nAC+FEL4ueX9cv7fBwF+q8/lxDfIOwP8mHS8qeS7ZXOBkADM7DNiDDAVAJapSu1TsAvx/NLFgZgPN7FNgAnB+1PUkmFkHYCAwOupaynGEmc01szfN7EdRF5OkG9DczN4zs1lmdnbUBZVlZo2AE/Af0FUW18WXrZznyl4N/QF4wMw+Aj4B5gBFNVxXVVSldimHmR2LB3nG2qIrE0J4GXjZzI4CbgP6RFxSwv3AtSGErWbl/ZWLzGx8PpBvzexE4BVg72hL+n/1gEOA3kBD4H0zmxlC+CzasrbRD5geQqjW/3riGuSLgN2TjjsCS5LfEEJYD5wHfoMR+KJki1qltcsPmdmBwFjgFyGEVVHXU1YIYZqZ7WVmrUIIcZh8KQ8YVxLirYATzawohPBKlEWV/LtM7L9hZg/H6DtbBKwMIWwENprZNOAgIE5BfjrVbFaB+DatfAjsbWZdzKw+/pt7LfkNZtas5DWAC4FpyX+JIvQacHZJ75XDgXUhhG+iLirOzKwT8BIwJE5XR2bWteQigZLeR/WBWPyQCSF0CSF0DiF0Bv4KXBp1iAOYWduk7+wwPGNi8Z0BrwK9zKxeSRPGT4AFEdf0/8ysKXA0Xme1xPKKPIRQZGbDgLfwO/P5IYR/mtnFJa+PBvYDnjSzrfid6AsyUZuZ/QW/M9/KzBYBNwM7JdX1Bt5z5XPgO0r+1xB1XWbWFigAdgWKzew3eE+gGv/hV4Xv7CagJfBwSQYUhQzMVFeFuk7BfyhvAb4HBiXd/Iy6tkhUoa5TgUvMrAj/zk6Py3cWQlhgZhOBj4FiYGwIYbvdOzNVW8nbBgJvl/yPoXqfn6HvWEREakhcm1ZERKSKFOQiIllOQS4ikuUU5CIiWU5BLiKS5RTkIiJZTkEuIpLl/g/6bl5sPA1DJwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(TT, LHS(TT), 'r-')\n",
    "plt.plot(TT, RHS(TT), 'b-');"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ahead-tribe",
   "metadata": {},
   "source": [
    "The intersection appears to be where it was exepcted.  The solution is confirmed!"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}